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The majority of current textbooks on device physics at the undergraduate level derive the diode equation 
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based on the diffusion of injected minority carriers. Generally the drift of the majority carriers, or the extent of 
drift, is not discussed and the importance of drift in the presence of a field in the neutral regions is almost totally 
ignored. The assumptions of zero field in the neutral regions and conduction by minority carrier diffusion lead 
to a number of pedagogical problems and paradoxes for the student. The purpose of this paper is to address 
the pedagogical problems and paradoxes apparent in the current treatment of conduction in the pn junction as it 
appears in the majority of texts. 
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I. INTRODUCTION AND BACKGROUND 



43 

The pn junction theory forms an integral part all physical electronics courses. At both undergraduate and graduate levels, the 
conventional analysis of the pn junction device under forward bias conditions follows closely Shockley's original treatment [1] 
in which the diode equation is derived based on the injection and diffusion of minority carriers. There are, however, a number 
of paradoxes in the treatment of the subject matter in the present textbooks (see below) which tends to mislead the students. 

q , Under zero and forward bias conditions, the pn junction displays a number of characteristic features with particular reference to 
carrier concentrations, exposed ionized dopants in the space charge layer (SCL) or the depletion layer, internal field, E(x) in the 
SCL, and the built-in voltage V . In a typical undergraduate course the arguments used in driving the diode equation follow a 

f~| sequence of simplifying assumptions: 

Qh| ( a ) The depletion region or the SCL (space charge layer) has a much higher resistance than the neutral regions so that the 
applied voltage drops across the depletion region. There is therefore no field in the neutral regions. 
' (b) The applied forward bias reduces the built-in potential Vo and allows the diffusion and hence injection of minority carriers. 
, From assumption (a), the law of the junction gives the injected minority carrier concentration in terms of the applied voltage. 
C<1 ' For example for holes injected into the n-side, 

in ; 

p n (0)= PnO exp(^) (1) 

where p n (0) is the hole concentration just outside the depletion region, at the origin x=0, p„o 1S the equilibrium hole concen- 
tration in the n-region, p n o=nj/Nd, V is the applied bias, and the other symbols have their usual meanings. 

(c) Since the electric field in the neutral regions is assumed to be zero, the continuity equation for the minority carriers is 
greatly simplified and becomes analytically tractable even at the junior undergraduate level. For holes in the n-region, under 
jj^. steady state conditions, 9p„/<9t=0 and the continuity equation is simply 
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-1^-^ = (2) 

e dx Thn 

• <-H , 

where J^„ is the hole current density and Thn is the minority carrier (hole) recombination time both in the n-region. It is tacitly 
assumed in eq.(|2j that the minority carrier concentration is much less than the equilibrium majority carrier concentration so that 
- - a constant minority carrier lifetime can be defined which is independent of the majority carrier concentration, n„. 

The hole current density however is simply the diffusion component as the electric field is assumed to be zero, 

Jhn = ~e-D hn —- (3) 
dx 

where Djj„ is the diffusion coefficient of the minority carriers (holes) in the n-region. 
Substituting eq. (0 into eq. (0, leads to 



(4) 
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and solving eq.0 for p„ one obtains, for a long diode, 



A Pn (x) = A Pn (0)exp(- — ) 

^hn 



(5) 



where Ap„=p n -p„o is the excess minority carrier concentration and x is measured from just outside the depletion region. The 
long diode assumption means that the length of the neutral region, 1„, is much greater than the minority carrier diffusion length, 
L/ m , defined as \/Dh n Th n - The long diode assumption however is not necessary for the derivation of the diode equation; it 
serves to simply the solution of eq.(@J to a single exponential rather than a hyperbolic sine function. 

With the minority carrier concentration given as in eq.(|5}, from eq.(|3} the hole current density is 



eD, m Ap n (0) x 
4n(x) = - exp(-- — ) 

There is obviously a similar minority carrier diffusion current density in the p-region for the injected electrons, i.e. 



(6) 



eL» ep An p (0) 



exp(- : 



(7) 
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where D ep and L ep are the electron diffusion coefficient and diffusion length, An p is the excess electrons concentration all in 
the p-region, and x' is distance measured away from the depletion region in the p-side. 

(d) It is assumed that the depletion region is so narrow that the currents do not vary across this region. Then the majority 
carrier current at x=0 is the same as the minority carrier current at x'=0. Similarly the majority carrier current at x'=0 is the same 
as the minority carrier current at x=0. Thus the total current density is 



J= J hn (0) + Jep(O) 

and using the law of the junction for Ap„(0) and An p (0) from eq^one obtains the diode equation, 



(8) 
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iD hn Ap n {0) e£> ep An p (0) 



Lhn 



(9) 



or 



J = e 



DhnPnO , D ep n p Q 
Lhn 



Jcp 



ex P (g) ^ 1 



(10) 



This is the general diode equation found in the majority of textbooks which follow the above sequence of steps, either through 
tacitly or explicitly stated assumptions in (a) to (d). Many texts, for simplicity, consider either the p + n or the n + p junction. For 
the p + n junction, Na >> No and eq.(l 1) becomes, 
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The assumption stated explicitly in (d) is commonly overlooked in many current textbooks on the subject and is one of the 
key assumptions in the derivation of equation Q as discussed, for example, by Moloney 1JJ] . 

Equation Q gives the impression that the conduction process in the p + n junction diode is the diffusion of injected holes in the 
n-region. The above steps invariably lead to a conclusion for the student that it is the minority carrier diffusion which constitutes 
the forward current. 

Close examination of the above steps in the derivation exposes a number of serious pedagogical paradoxes and problems for 
the student and the instructor. The diode equation in eq. [^is so entrenched in our teaching of the pn junction that it has been 
used to design many simple but fruitful laboratory experiments as reported in various journals. 
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II. PEDAGOGICAL PROBLEMS 



The conventional undergraduate level treatment in Section I leads to a number of pedagogical problems and paradoxes. We 
cite those we have encountered frequently in two undergraduate classes during the treatment of the forward biased long diode: 

(a) If there is no field in the neutral region then there can be no net charge at any point in this region inasmuch as 
dE/dx = 0. Then the excess majority carrier concentration should follow the decay of the excess minority carrier concen- 
tration, An„(x)=Ap„(x), which then follows eq. [5] But the gradients of An„(x)=Ap„(x) along x must be the same as well. 
Therefore majority carriers diffuse towards the right as well and since in silicon D e > Dh, and electrons are negatively charged 
the net current is actually in the reverse direction. The current must be in the opposite direction to the applied voltage! 

(b) The minority carrier current, which is due to diffusion, decays with x. Since the total current must be constant, the majority 
carrier current must increase with x. However, there is no field in the neutral region which means that electrons must diffuse. 
From the first paradox above this diffusion can not make up for the decay in the hole current. 

(c) Far away from the depletion region, both the hole and electron concentrations are almost uniform. If there is no electric 
field, then the current, due to diffusion, must vanish. How is it that the current stays constant in the neutral region (indeed in the 
whole device)? 

(d) The absence of an electric field in the neutral region means that An„(x)=Ap„(x). But when holes are injected into the 
n-region they recombine with electrons so that , intuitively, the majority carrier concentration should decrease not increase. 

The concept of field free neutral regions is so deeply rooted in the present treatment that many authors do indeed show the 
excess majority carrier concentration increasing towards the junction like the excess minority carrier concentration. This fact 
alone is contrary to student intuition that, if anything, the excess majority concentration should remain uniform. It is clear that 
the instructor has a responsibility to clear this paradox. It is interesting to note that a large number of authors sketch the carrier 
concentration profiles with the majority carrier concentration shown as uniform whereas others show the excess majority carrier 
profile following the excess minority carrier concentration profile as indicated in Table I for a survey of a large number of books 
on the subject. The differences in the diagrams only add confusion to the student's understanding. Few authors allude to the 
presence of the electric field and majority carrier drift to overcome some of the problems listed above. There is however no 
satisfactory treatment in the majority of the textbooks used in English speaking countries we have examined as illustrated in 
Table I. In most books the field in the neutral region is totally neglected in the treatment. The result, we believe, is pedagogical 
paradoxes and a student who is confused. It seems that only some of the early books published in the sixties consider the need 
and the effects of the field in the neutral regions in their discussion of conduction in the forward biased pn junction. 



For simplicity we consider a long p + n junction diode. Equation (|2ji describes its conventional current-voltage characteristics. 
We use the parameters listed in Table 2 which represent "typical" parameters for a p + n junction Si diode albeit classroom 
parameters. We also assume small injection so that p n o << n n Q (or Nd). The latter assumption means that the minority carrier 
recombination time remains a useful parameter . For forward bias we take the voltage across the diode to be typically 0.55 V. The 
depletion region extends essentially into the n-side and its width W, is much shorter than the hole diffusion length in this region 
as indicated in Table II. Similarly the relatively tiny extension of the depletion width into the p + region is much shorter than the 
electron diffusion length there. Lengths of the neutral regions are taken to be about ten times their minority carrier diffusion 
lengths to represent a "long-diode", i.e. l n =10L/, and l p =10L e . 



The first attempt to overcome the problems listed in Section II is to allow some of the applied voltage, a small fraction of it, 
to drop across the neutral n-region of the p + -n junction. This is easily accepted by the student since the neutral regions must 
have some finite resistance even though much smaller than the depletion region. This means that the law of the junction remains 
approximately valid. What is the field in the n-region? 

The total current through the p + n junction must be continuous. This means that at any point in the n-region, 



III. DIFFUSION OR DRIFT? 
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where E„ is the field in the n-region at x. Initially E n is assumed to be small but finite. 

Since n n >> p n , n n w No (small injection), and An„(x)=Ap„(x) the above equation simplifies to, 
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J = e(D en - D hn ) 



dAp v 
dx 



ep ien N D E n (x) 



(14) 



The requirement of an internal field is quite transparent from eq. The minority concentration gradient is negative and 
D en > Dhn which means that the first term in eq.(|3} is negative so that current is in the negative direction. Unless there is an 
internal field drifting the majority carriers it is not possible to obtain a positive current. 

The pedagogic development at this point must make use of the excess minority carrier concentration in eq|5] If the field is 
indeed sufficiently small it may be assumed that the excess minority carrier profile, Ap„(x), is still given by eq|5] The validity 
of this assumption will be demonstrated below with an illustrative example. One can also assume that eq. (|2j can still be used to 
describe the total diode current. Then from eqs. |3]|2]and|3] one can obtain the field E„ in the n-region 
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kT 

&Lh.r) 



1 eV 



l + (6„-l)exp(— — ) 
J-'hn 



(15) 



where we have used the definition b„ = /i en //i/j n =D e „/D/j„(> 1). Equation (0} describes the field outside the SCL in the 
so-called neutral region that is needed to maintain the diode current. In the n-region the field increases towards the SCL. The 
increase in the field is required to make up for the negative electron diffusion current. Far away from the depletion region the 
current is maintained by a constant field of magnitude, 



kT 

6-Lhn 



1 eV 



(b n - l)exp(-- — ) 

J^hn 



(16) 



An interesting feature is that the magnitude of the field increases exponentially with the applied voltage contrary to student 
intuition based on the applied voltage simply dividing between the resistance of the depletion region and the resistance of the 
neutral region. 

With the field given in eq[0]a paradox mentioned in Section II develops in that E varies spatially across the neutral region so 
that dE/dx is not zero. Gauss equation in point form (or the Poisson equation) in the n-region is: 



E n (x i — * oo) = 



kT \ l_ eV 
\N D J \eL hn ) b n eXp{ kT } 



(17) 



f^L = A p n (x) - An n (x) (18) 
e dx 



where e is the total permittivity of Si (=e e r ). 

Since Ap„(x) is determined by eq|5] the excess majority carrier concentration is: 



An n (x) = Ap n (x) 



e dE n 
e dx 

Ap n (x) 



Ap n (x) 



e(b n - l)kT 
e 2 b n N D L\ n 



(19) 



Substituting typical values for e, b„, No, L/ m from Table II into eq.(|2j shows that the second term is ^4.1xl0 -6 . Thus 
An„(x)=Ap„(x) and the charge neutrality condition for all practical purposes remains valid. We have found the requirement of 
An„(x)=Ap„(x) to be somewhat contrary to student intuition. This is further exasperated by many texts showing the majority 
carrier concentration uniform in the neutral regions (see Table I) which misleads the student. Qualitatively, the injected holes 
into the n-region disturb the charge neutrality and set-up a field here which then drives the electrons towards the SCL until a 
steady state is reached between electron drift and diffusion. Thus the increase in the majority carrier concentration towards the 
SCL is due to the driving effect of the field, E„, even though it appears at first that An„ should decrease towards theSCL as 
injected holes recombine with electrons 

Once the field in the n-region is given as eq^l the student can readily calculate the various contributions to the total current 
density using eq^] The magnitudes of the various current components (majority carrier diffusion and drift, and minority carrier 
diffusion and drift) and their directions are listed in Table II. In general, the drift of the majority carriers is the most significant 
contribution to the pn junction diode current. How is then the diffusion terminology comes to appear in explaining the diode 
current even though the biggest contribution is drift? 
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Given that the depletion region is very narrow and that recombination in this region is negligibly small due to the very small 
concentrations of carriers, then in the steady state one must have = and the electron current J e must be constant through 
the SCL. But, electron drift at x'=0, i.e minority carrier drift, is negligible and the electron current there is primarily a diffusion 
current just like the hole current in the n-region. Thus the total electron current at x=0 must equal to the electron diffusion 
current at x'=0. Similarly the hole diffusion current at x=0 is equal to the total hole current at x'=0 which is essentially by drift. 
It is apparent that by evaluating the minority carrier diffusion current just outside depletion region we are indirectly determining 
the total majority carrier current at the other side of the depletion region. This is a subtle point that seems to have been short 
circuited in the majority of the texts. Consequently, the diode equation stated in equation (10) is only valid if the SCL width is 
much shorter than the minority carrier diffusion length in that region 

An important paradox that must be addressed by the instructor is the much cherished minority carrier concentration profile 
stated in eq|5]for a long diode. Equation|5]is the solution of the continuity equation in the absence of an electric field. As a first 
step one can assume a constant field, E n , in the neutral region to examine its effects on the excess minority carrier concentration 
profile. In the presence of a constant field, the general continuity equation in eq|2]leads to 



Uhn o H , hn-*- J Ti 
QX 



dPn _ Pn_ 
dx Thn 







(20) 



Since this is a linear differential equation, the undergraduate student can readily solve it or accept its solution by substituting 
the solution into eq.(|3}- For a long neutral region, the solution is: 



Ap n (x) = Ap n (0)exp( ) 



(21) 



where a„ is a "length constant" which can be determined by substituting eq.@ into 0; 



(Thn^hnE n )a n - {D hn T hn ) = 



(22) 



Solving this quadratic equation we find, 



\J (DhnThn) [Va n + 1 + Ctn] = L hn [\/ a n + 1 + . 



(23) 



where a„ is defined as 



drift 



Schubweg 



1 



2 Dif fusion Length 



(24) 



The parameter a n represents the comparative effect of drift to diffusion since L<j r j f t = 1 liThnl l 'hr^n is the so-called Schubweg 
of the minority carriers, distance drifted before recombination. If the field is small a „will be small and in the limit of zero 
field, E„ ~ 0, a n ~ Lhn an d the theory approaches the conventional zero-field treatment. At a forward bias of 0.55V, the field 
is maximum at x=0, and using this maximum value one finds a n =0. 00106 and a„=(1.001 l)L/ m . At V=0.65V, a„=(1.05)L/,„ 
and a„ is still very close to L/ m (within 5%) even though the injected hole concentration is now only ten times smaller than the 
equilibrium majority concentration which sets the limit of small injection. Although the solution in eq.(|3} does not apply when 
the field is non-uniform as in eq[21 it does nonetheless provide convenient means for the student to examine the possible effect 
of the field on the excess minority carrier profile. 

The field in the p-region can be similarly derived. The total current in the p + region is 



J = eD e 



dAn r 
' dx 



dApp 
dx 



euhpPpEp = constant 



(25) 



Since the total current must be constant and assigned to be described by eq. ^] using the corresponding version of eq. |5]for 
minority and majority carrier excess concentration in the p-region one can derive 



eV 

E p (x) = cxp( — ) 



kT 



kT' \ N D N A \eL h 



/ Hhn \ 



1 - C p exp(- — ) 



(26) 
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where C p is defined as 

in which b p = /J, ep /[ihp is tne electron to hole drift mobility ratio in the p-region. We assumed that, as usual under forward 
bias, V >> kT/e. Equation (|9jl shows that the field is minimum right at the SCL, x'=0, and increases exponentially to a 
constant value away from the junction. Substituting typical values from Table II shows that E p « E n . When the two fields 
are compared, one finds that E n is at least three orders of magntiude greater than E p . In fact E p is almost unifrom in the p + - 
region. Most interestingly and importantly, even though E p is even smaller than E„, its effect is most significant. One readily 
can calculate the contributions of each term to the current density in the p-region from eq. [24] The values at the SCL are listed 
in Table II where it is apparent that the current is carried almost totally by majority carrier drift. A distinctly different behavior 
in the p + -region from that observed in the n-region is the fact that the majority carrier diffusion is insignificant and that minority 
carrier diffusion, though larger than majority diffusion, is some three orders of magnitude smaller than majority carrier drift. 
From the above discussion for conduction on the n-side it is clear that in deriving the diode equation we are calculating the hole 
(majority) drift current in the p + -region by evaluating the hole (minority) diffusion current in the n-region simply because the 
total hole current does not change through the SCL as long as the latter is thinner than the hole diffusion length. 

It is always useful for the student to reconfirm that the majority of the voltage drops across the depletion region by evaluating 
the voltage drop across the neutral regions. If V„ is the voltage drop across the n-region then 

V„ = / " E n dx (28) 

or 

/ //, y ( kT 



Vn = vt) VuTJk ^ + ^ - ^ exp ^ (29) 

Equation l|2} shows that the voltage drops increases exponentially with the applied bias contrary to an intuitive guess. Using 
typical values, at V=0.55V, V„ is 0.00168V, whereas at V=0.6V, V„ is 0.019V and the injected hole concentration in the n-region 
is about 1 1% of No which is the limit of small injection. At V=0.65V, V n becomes 0.121V which is quite significant but at this 
bias voltage the injected hole concentration is no longer small compared with No. There is a clear indication that as the voltage 
across the diode increases more and more of the applied voltage drops across the neutral regions which deteriorates the law of 
the junction. It is not difficult to show that since E p << E n and l p « l n , the voltage drop across the p + -region is orders of 
magnitude smaller than V n . 

IV. NUMERICAL DESCRIPTION OF THE PN JUNCTION: 

While the previous sections discussed the Diffusion/Drift approximation problem, this section is oriented towards an exact 
description of the PN junction without making any assumptions. 



Starting from the constitutive system of equations 130]: 



dib 

~T = -E (30) 

dx 



dn di/j 

J n = e[D n - fJ. n n— \ (31) 

dx dx 



d P di; 

J o = - e[D *te +flpP te ] (32) 

1^ =R(n,p) + G(x) (33) 
e dx 
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-^ = -R(n,p)-G(x) (34) 
e dx 



with the recombination term of the Shockley-Read-Hall form: 



R(n>p) = ("P-"?) (35 ) 
V ,y) T in + T 2P + T 3 



where Ti, T2and T3are time constants. Typically for Silicon Ti = t p = 10~ 5 sec,T2 = r„= 10~ 5 secand T3= (r p + T n )r\i 13111 . 

In order to solve the boundary value problem associated with the above system (when a voltage is applied to the PN junction) 
we transform it into a new hybrid system of first-order (Current and carrier density equations) and one second-order differential 
equation (Poisson equation). 

The mathematical/numerical reasons for performing this transformation reside in the fact the above system is a "singularly 
singular perturbed problem" |38l |39[. Many algorithms |3SJ |4l[ have been developed in order to deal with this difficulty 
stemming from several facts: 



1. ip, n and p are fast variables in comparison with E, J„ and i p 13 W32U33H . 

2. Near the limits of the depletion layer the values of n and p change by several orders of magnitude making the space-charge 
zone a double boundary layer. This difficulty is of the same type as the one encountered in Hydro/Aerodynamics where 
the fluid velocity field changes by several orders of magnitude near an obstacle. 

Recognizing the difficulty due to the presence of the space-charge layer, a standard way to find a valid solution is to treat 
the boundary layer separately from the rest of the diode. In spite of the success of this approach I36ll37ll . one might feel 
uneasy about this artificial dichotomy and rather tackle the problem with new powerful mathematical/numerical methods 
that handle the layer and the rest of the device on the same equal footing. 

3. When the above system is rewritten explicitly in terms of the Poisson equation as we will do below, the second spatial 
derivative of the electric potential is multiplied by a very small number A 2 (A ~ 10 -4 tol0 -3 ). Actually, this is the reason 
the problem is called singularly perturbed: the solution with A = is entirely different from the solution with A finite but 
small il. 

One of the early algorithms aimed at circumventing the above difficulties is the Scharfetter-Gummel algorithm 1 34]. The latter 
attempts at segregating the fast/slow variables by integrating out the fast variables over some small interval while holding the 
slow variables constant over that same interval. The Scharfetter-Gummel algorithm leads to a spatial exponential discretization 
that will alleviate for the rapid variation of the fast variables. 

Many variants of the Scharfetter-Gummel algorithm have been developed [35] in order to cure some of its shortcomings 
which generally are numerical oscillations and crosswind effects. These lead to a loss of accuracy of the solution and sometimes 
preclude convergence towards the solution. 

We decided not to use the Gummel algorithm or any of its variants but rather to tackle the problem from the singular perturba- 
tion point of view since this approach is more rigorous and leads to a better understanding and control of the instability problem 
encountered in the semiconductor system of equations. 

We first transform the system in the following dimensionless two-point boundary value problem with no generation processes 
considered: 

dn dtp 

— = Ci J n + n— 
dx dx 

dp dip 

— = -6 2 J P -p— 
dx dx 

dJ„ (np - 1) 

— ^3" 



dx n + Tip + r 2 



dJ p (np - 

— — O; 



dx n + Tip + t 2 
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= C 4 (n-p + N D ) 



The constants Ci,C2,C3 and C4 are given by: 

_ JqLd r , _ JqL d _ eL D n t _ eL 2 D m 
Gl = FT' ° 2 — FT' ° 3 — t rri — — p7 — 



where the Debye length Lp is given by — y k ^J s and the scaling current J = "' m £^ bT . 

The thermal voltage UT=ksT/e (T is the temperature and ks is Boltzmann constant). The time constants are now Ti=T2/Ti 
and T 2 =T 3 /(T2ni). 

The above system is now in the appropriate form to integrate with a powerful B-spline collocation based algorithm specifically 
tailored for two point singularly perturbed boundary value problems: COLS YS I39tl40l EH. The algorithm is based on a control- 
lable meshing technique I4ll 14211 of the boundary layer which will lead essentially to the damping of the existing singularities. 
The layer-damping mesh, being exponential in nature, encompasses the Scharfetter-Gummel case and can be shown rigorously 
to have the form: 

hi = hi^i exp(a/ii_i/A) 

where hi is the i-th mesh point, a is a constant related to the required accuracy and the nature of the collocation and Ags the 
singular value parameter. 

Previously, Markowich et al. |38] tackled this problem from the same angle but they solved the symmetric diode case with 
one boundary-layer at the origin. 

In this work we tackle the non-symmetric case where a double-boundary layer is present around the origin starting from very 
accurate initial conditions. 

Varying the applied bias by steps of 0.1 V we calculate the carrier density profiles, the potential and the electric field. 

The drift and diffusion current density profiles are also obtained. Several tests are used in order to check the validity of the 
solution obtained. The first test is an accuracy test whereby we require a given accuracy and check whether the criterion is met. 
The next test is based on the requirement of convergence: the collocation builds a non-linear set of equations that has to be 
solved iteratively. The additional tests are the independant checks of the constancy of the current densities locally on each side 
of the junction and globally over the entire junction. The tests are shown in the current density profile figures. 

The final test we use is the approximate validity of Shockley's equation. Varying the voltage, we obtain the IV characteristics 
of the junction and we compare it to the Shockley's case. Since we have used the Shockley-Read-Hall recombination term all 
over the PN junction we do not expect the Shockley's case but rather the general form: 

/ = Is[exp(eV/r)k B T) - 1] 

The comparison of the obtained IV characteristic to the Shockley formula is displayed in Fig. 1. The calculated characteristic 
falls between the two Shockley curves J7 = 1.1 and i] = 1.2 in a finite current interval. This means, the general Shockley formula 
is not valid, within the singular perturbation annroach for arbitrarv current values 




FIG. 1: IV characteristics obtained for ther PN junction and Shockley approximation with 77 = 1.1 and 77 = 1.2. 
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V. CONCLUSIONS 

The physics of the PN junction is gaining back interest with new developments in the area of nanoelectronics specially in the 
area of spintronics where one has to account for the spin of the carriers in addition to their charge. The usual approximations that 
are valid and successful in the description of the PN junction physics at the micron scale must be entirely reviewed and adapted 
to the nano scale. The diffusion/drift approximation as well as the nature of the singularities of the problem have been reviewed 
and reformulated in a way such that the underlying assumptions are revealed with their consequences. 
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VI. TABLES 



AUTHOR 


REF. 


SEMIQUANTITIVE 
DRIFT ANALYSIS 


FIELD IN NEUTRAL 
REGIONS 


MAJORITY CARRIER 
CONCENTRATION 


COMMENT 


Zambuto [2611 


Ch. 6 


Qualitative 


Not discussed 


Not shown 


Undergraduate 


Yang |21] 


Ch. 4 


Qualitative 


Not discussed 


An n (x)=Ap n (x) 


Undergraduate 


Elliott and Gibson |5] 


Ch. 9 


None 


Not discussed 


n n shown constant 


Undergraduate 


Fraser (UK)|7] 


Ch. 3 


None 


Not discussed 


n„ shown constant 


Undergraduate 


Seymour 1 1711 


Ch. 3 


None 


Not discussed 


Not shown 


Undergraduate 


Streetman l20ll 


Ch. 5 


Some (Example 5.4) 


Mentioned but not dis- 
cussed 


Not shown 


Undergraduate 


Valdes [22] 


Ch. 9 


Yes 


Discussed 


An„(x)=Ap„ (x) 


Out of print / Undergraduate 


Sze | 27] 


Ch. 3 


None 


Not discussed 


Not clear 


Undergraduate 


Carroll (UK) H 


Ch. 4 


None 


Mentioned but not dis- 
cussed 


Not clear 


Undergraduate 


Solymar and Walsh 
(UK) 1 19] 


Ch.9 


None 


Not discussed 


An n (x)=Ap n (x) implied 


Undergraduate 


Shur 1 18] 


Ch. 2 


Some (Fig. 2.3-8) 


Not discussed 


An n (x)=Ap n (x) implied 


Senior UG / Graduate level 


Pulfrey and Tarr 
(Canada) 0] 


Ch. 6 


None 


Not discussed 


An n (x)=Ap n (x) implied but 
n n =constant in diagrams 


Undergraduate 


Colclaser and Diehl- 
Nagle H 


Ch. 7 


None 


Not discussed 


Not mentioned and not shown 


Undergraduate 


Navon 1 14] 


Ch. 6 


None 


Not discussed 


An n (x)=Ap n (x) in diagram 


Undergraduate 


Sze | Si 


Ch. 2 


None 


Not discussed 


Not clear 


Senior Undergraduate / Grad- 
uate 


Wang 1 25] 


Ch. 14 


None 


Not discussed 


Inferred 


Senior Undergraduate 


Goodge (UK) |8] 


Ch.l 


None 


Not discussed 


An n (x)=Ap n (x) in diagram 


Undergraduate 


Grove 1 12] 


Ch.3 


None 


Not discussed 


Not clear 


Out of print / Undergraduate 


VanDerZiel 1 23] 


Ch. 15 


None 


Not discussed 


Not clear 


Out of print / Undergraduate 


Lonngren I13ll 


Ch. 6 


None 


Not discussed 


Not discussed 


Undergraduate 


Tyagi Ei] 


Ch. 7 


Qualitative 


Mentioned but not dis- 
cussed 


An„(x)=Ap„(x) (Fig.7.3) 


Undergraduate 


Gibbons \9] 


Ch. 6 


Yes 


Discussed 


An n (x)=Ap n (x) 


Out of print / Undergraduate 


Ferendeci |6J] 


Ch. 8 


None 


None 


Not clear 


Undergraduate 


Allison |2] 


Ch. 7 


None 


None 


Not clear 


Undergraduate 


Neamen |ll_5J] 


Ch. 8 


Some 


E estimated 


Not clear 


Undergraduate 


F.Wang 1 24] 


Ch. 8 


None 


None 


Not clear 




Gray et al. Il ill 


Ch.2 and 
App. B. 


Some 


Yes 


An, I (x)=Ap, I (x) in diagram 


Out of print / Undergraduate 
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PROPERTY / PARAMETER 


TYPICAL VALUE 


COMMENT 


Permittivity 


11.9 




Intrinsic concentration, ni, cm -3 


1.5xl0 10 




Donor concentration, No, cm -3 


5xl0 15 




Acceptor concentration, Na, cm -3 


10 19 


p + -n junction 


Equilibrium hole concentration 
in n-region, p n0 , cm" 3 . 


4.50 xlO 4 




Equilibrium electron concentration 
in p-region n p o cm 3 


22.5 




Hole recombination time in n-region, Thn, S 


5x10" ' 
hn 1 + 2X 10- 17 JV D 

Thn = 4.54 x 10~ 7 


T en = Thn 


Electron recombination time in p-region, r eP; s 


5xl0~' 
ep l + 2xlO- 17 W^ 

T cp = 2.49 x 10~ 9 




Electron drift mobility in the n-region, 
cm 2 V 1 s 1 . 


Hen - 88 + 1252 18 

^ l + 6.984xlO- ls N D 
Hen = 1298 


b ~ = m; = 2 - 86 


Hole drift mobility in n-region, 
„ m 2 v -i -l 
cm V s 


a hn - 54.3 + 407 „ 

^ nn 1 + 3.745X10- 18 JV D 

ft h „ = 453.8 




Electron drift mobility in the p-region, 
cm 2 V -1 s -1 . 


Mep 88 + 1 + 6 . 984x10 -18 JVa 
llep = 105.7 


b " = ^ = 163 


Hole drift mobility in p-region, 
cm 2 V~ 1 s~ 1 . 


Mdp - 54.3 + 1+3 745x10 -is Na 
(J-hp = 64.9 




Electron diffusion coefficient in n-region, D en , cm 2 .s~ 1 


33.55 


Den _ kT 
fJ-pn e 


Hole diffusion coefficient in n-region, D hn , 
cm 2 .s -1 


11.73 




Electron diffusion coefficient in p-region, D ep , cm 2 .s _1 


2.73 




Hole diffusion coefficient in p-region, D^ p , cm 2 .s -1 


1.68 




Electron diffusion length in n-region, L en , cm 


3.90x10- 3 




Hole diffusion length in n-region, Lhn, cm 


2.31 xl0~ 3 




Electron diffusion length in p-region, L ep , cm 


8.24x10- 5 




Hole diffusion length in p-region, L^p, cm 


6.46x10- 5 




Builtin potential, \ b i, V 


0.854 




E bi , V.cm -1 


3.60 xlO 4 


No bias 


Width, W, of depletion region, cm 


4.741 xl0~ 6 


On n-side. Much shorter than hole diffusion lengt 


Width, W n , of depletion region in n-side, cm 


4.739xl0~ 5 


Much shorter than hole diffusion length in n-regio 


Width, W p of depletion region in p-side, cm 


2.369X10- 8 


Much shorter than electron diffusion length in p-r 


Length of n-side 


2.31 xlO" 2 


10Lh Long diode 


Length of p + -side 


8.24xl0~ 4 


10L e Long diode 


FORWARD BIAS, V 


0.55 




Built-in electric field, V.cm -1 


2.15X10 4 


Smaller than zero bias case 


Width, W, of depletion region at 0.55V, cm 


2.83xl0~ 5 


On n-side. Narrower under forward bias. Much 
side. 


p n (0), injected hole concentration at x=0 


7.81 xlO 13 




rip(0), injected electron concentration at x'=0 


3.91 xlO 10 




Pn(0)/N D 


0.0156 


1.56%, small injection. At V=0.60V, this become* 


n p (())/N A 


3.91 xlO" 9 


Extremely small injection 


J 0h , A.cm -2 


3.663 xlO" 11 




J 0e , A.cm -2 


1.195X10" 13 


Joe << Joh 


J , A.cm -2 


3.675X 10" 13 


~Joh 


Jo. 55, A.cm -2 


0.0638 


.7 = .7 [oxp(y/fcT) - 1] 


E nx , V.cm -1 


0.0613 


Field in n-region far away from junction 


Enmax, V.Cm _1 


0.1749 


Field just outside SCL 


E nma x/E bi 


8.13X10" 6 




V^at V=0.55V 


0.00168 


Very small, V„is 0.3% of bias 


V^at V=0.60V 


0.0116 


V„is 12.3% of bias 


Epmax = E pX , V.cm -1 


6.1 xlO" 4 


Field far away from junction 


Epmax^bi 


2.8xl0" 8 


Extremely small 


v P 


5.04X10" 7 


Extremely small, V p < < V n < < V 


JUST OUTSIDE SCL ON n-SIDE, x=0 






Majority drift current, A.cm -2 


0.1820 


Largest magnitude in positive direction 


Majority diffusion current, A.cm" 2 


-0.1818 


Opposite direction, about the same magnitude as r 


Minority drift current, A.cm " 2 


0.001 


Smallest magnitude 


Minority diffusion current, A.cm -2 


0.06359 


About a third of the magnitude of majority drift. 


JUST OUTSIDE SCL ON p-SIDE, x'=0 






Majority drift current, A.cm -2 


0.06351 


Largest magnitude. Dominates conduction. 


Majority diffusion current, A.cm - 2 


-0.000891 


Opposite direction 


Minority drift current, A.cm " 2 


2.33X10" 19 


Smallest magnitude-virtually zero 


Minority diffusion current, A.cm -2 


0.00255 


Next largest magnitude 



